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Abstract 

Based on the renormalisability of the SU(n) theory with massive gauge bosons, we start with the 
path integral of the generating functional for the renormalized Green functions and develop a method 
to construct the scattering matrix so that the unitarity is evident. By using as basical variables the 
renormalized field functions and defining the unperturbed Hamiltonian operator Ho that, under the 
Lorentz condition, describes the free particles of the initial and final states in scattering processes, we 
form an operator description with which the renormalized Green functions can be expressed as the 
vacuum expectations of the time ordered products of the Heisenberg operators of the renormalized 
field functions, that satisfy the usual equal time commutation or anticommutation rules. From 
such an operator description we find a total Hamiltonian H that determine the time evolution of 
the Heisenberg operators of the renormalized field functions. The scattering matrix is nothing but 
the matrix of the operator U(oo, — oo), which describes the time evolution from -co to oo in the 
interaction picture specified by H and Ho, respect to a base formed by the physical eigen states of 
Ho- We also explain the asymptotic field viewpoint of constructing the scattering matrix within our 
operator description. Moreover, we find a formular to express the scattering matrix elements in terms 
of the truncated renormalized Green functions. 

PACS numbers: 03.65.Db, 03.80.+r, 11.20.Dj 
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I. Introduction 



Although the SU(n) theory with a mass term of the gauge fields was considered to be in general 
nonrenormalisable for a long time, it was repeatedly studied (see for example Refs. [1-8]). In particular, 
we have proven the renormalisability of the theory under the original expression of the mass term of the 
gauge fields [8]. Therefore it becomes possible to to deeply studied the scattering matrix of such theories. 
In the present paper, instead of checking the negative arguments concerning the unitarity we will develop 
a method to construct the scattering matrix with the renormalized field functions as basical variables so 
that the unitarity is evident. 

Owing to the renormalisability of theory, one can start with the path integral of the generating 
functional for the renormalized Green functions and form an operator description with the renormalized 
field functions as basical variables. This means that the Hciscnbcrg operators of the renormalized field 
functions satisfy the usual equal time commutation or anticommutation rules and the renormalized Green 
functions can be expressed as the vacuum expectations of the time ordered products of such Heisenbcrg 
field operators. On the other hand starting from the path integral of the generating functional for the 
unrenormalized Green functions, one can form another operator description with the bare field functions 
as basical variables. It should be emphasized that such two operator descriptions are of different types. In 
particular, the field operators belonging to the two operator descriptions operate on different state spaces 
and can not be transformed into each other by unitary transformations or by scale transformations. 

According to the very idea of renormalisation we decide to choose the first kind of operator de- 
scription. With the renormalized field functions and the renormalized parameters we will define a non- 
interacting Lagrangian C which, under the Lorentz condition, describes the free particles of the initial 
and final states in scattering processes. Next we will find an unperturbed Hamiltonian H operation in a 
way used in the Gupta-Bleulcr quantization and define a Quasi-Gupta-Bleuler subspace (QGB) which is 
spanned by all the physical eigen states of H . By giving an operator expression for the generating func- 
tional for the Green functions of the free fields, we will able to do the same for the generating functional 
for the renormalized full Green functions. Then we can express the renormalized Green functions as the 
vacuum expectations of the time odered products and find the total Hamiltonian H that determine the 
time evolution of the Hciscnbcrg operators of the renormalized field functions. With the QGB subspace 
and the operator U(oo, —oo), which describes the time evolution from — oo to oo in the interaction picture 
specified by H and H , we can easily define the scattering matrix and make the unitarity evident. 

The operator description will be given in section 2. In section 3 we will explain the definition and 
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the unitarity of the scattering matrix and briefly describe the asymptotic field viewpoint based on our 
operator description. In section 4 a formular will be found to express the scattering matrix elements 
in terms of the truncated renormalized Green functions. Concluding remarks will be given in the final 
section. 

II. Operator description with renormalized field functions as basical variables 

Assume that the theory has been renormalized under a proper renormalization scheme so that the 
renormalized parameters such as the masses represent the physical ones. The generating functional for 
the renormalized Green functions can be expressed as 

2^[J,x,x,V,v] = ^J II 6(d»A a ,„(x'))e X p{i(l eS + I s )}, (2.1) 

a' ,x' 

where A afl ,uj a ,ZJ a and tp,tjj stand for the renormalized gauge fields, ghost fields and material fields, all 
the parameters are also renormalized ones, N is a constant to make [0, 0, 0, 0, 0] equal to 1, I s is the 
source term 

Is = J d i x[j^(x)A all (x) +X a (x)Ua(x) + UJ a (x)Xa(x) 
+Va( x )'<Pa(x) +1p a {x)T] a (x)] . 

The action I e g contains the counterterm and is formed by the following Lagrangian 

2eff = ~F aiM ,F^ + X -M 2 A a ^ a + + C^ A 

+ ( - d ll U a {x))D^ b U3 b {x) + Ccount ■ (2.2) 

The free evolution of the physical particles of the initial and final states in scattering processes should be 
described by combining the Lorentz condition with the Lagrangian 

£>(*) = -\4%F^ + \M 2 A a ,{xW a (x) + C^x) + ( - d^ a (x))D^ b (x) , (2.3) 

where 

= d A - d A 

In order to describe such free particles by an operator method, we replace (2.3) by 

Cq(x) = C GB (x) + ^M 2 A a( t(x)AZ(x) + C*(x) + ( - d^ a (x))D^ b {x) , (2.4) 
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where £gb(x) is a Gupta-Bleuler Lagrangian 

C GB (x) = -\Fi%F^- 1 -(d^A a ,) 2 . 

Then we perform a canonical quantization with Co(x) and define the corresponding Hamiltonian operator 
H and field operators. Moreover we have to choose the so-called physical states and define a physical 
initial-final state subspace according to Lorentz condition and the residual symmetry. This subspace 
is spanned by all the physical eigen states of H and can be formed in a similar method used in the 
Gupta-Bleuler quantization although there are differences caused by the gauge field mass term (see the 
Appendix or Ref. [9]). We will simply call such a subspace as Quasi-Gupta-Bleuler or QGB subspace. 
If we can find a total Hamiltonian which determines the time evolution of the Heisenberg operators of 
the renormalized field functions, we will be able to construct the scattering matrix with the help of the 
QGB subspace and the time evolution operators [/(oo, — oo) which describes the time evolution from — oo 
to +00 in the interaction picture defined by the total Hamiltonian and Hq. It will be seen that these 
tentative ideas are correct. However, we will have to introduce the £— gauge fixing term and construct 
an appropriate operator expression for the generating functional for the renormalized Green functions. 
With the help of the the gauge fixing term Cgf = — 5f (^ M Aif*) one can express (2.1) as 

Z [R] [J, % X, V, v] = lim Zf ] [J, x, X, V, v] , (2.5) 

where 

z\ m [J,X,X,V,v] = ]^/ ^,w,o;, ; 0,V]exp{i(/ eff + 7 s )}, (2.6) 

where N$ is a constant to make Z^ [0, 0, 0, 0, 0] equal to 1, and 

Jeff = Jeff + J d 4 xC GF (x). (2.7) 

Needles to say, quantities defined by using Zp' [ J, x, X, V> v] depend on £, but this parameter will often 
be omitted. 

Let Z [J,x,X:Vi 7 l] be the generating functional for the Green functions defined by the effective 
Lagrangian Co, namely 

Z°[J,x,x,V,v} = ^- J T>[A,LJ,u>rf,il>] J] S(d»A a> (x'))cxp[i( J d 4 x£ (x) + 1.) } , 

a' ,x' 

where 7V is a constant to make Z o [0, 0, 0, 0, 0] equal to 1. One thus has 
Z°[J,x,X,V,v] = (0|Texp{iy d i x[J^(x)A a „(x) 

+ X a (x)u a (x) +uj a {x)Xa{x) +rj a {x)^ a {x) +lp a {x)r) a (x)]}\0) , (2.8) 
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where A atl {x),Lu a (x),ip a {x),- ■ ■ are field operators satisfying the usual equal time commutation or anti- 
commutation rules and their time evolution are determined by H . The state |0) is the eigen state of H 
belonging to the lowest eigenvalue (assigned to be zero). Define Hi so that 

J eff = - J d i xH I (x) + J d 4 xC (x) . (2.9) 

Then one gets from (2.6) 

f • f ,4 jr. ( 5 5 5 6 5 \ 1 z r T _ _ i 

ocexp<^ -i / d xHi — , , . u , 7j=-, r , 7 ., x , . >Z [J,x,X,V,V] 

{ J \iSJ(x) iS x (x) (-i)Sx(x) i6rj(x) (-i)Sr)(x) J J 

cx(0|T(cx P {i J d 4 x[JZ{x)A afl (x)+x a {x)u a (x) 

+LJ a (x)xa(x)+rj a (x)Mx)+^ a (x)Va(x)]}Te^-i J AH,W))|0). (2.10) 

where Hi ( ) is a differential operator obtained from Wj(a:) by replac- 

ing the field functions A aii (x), w a (x), uj a {x), ip aa {x) and i/) a p(x) with iSJ i (x) , wijx)' =I5^R ' ^l(^) 
and zig^~(^) respectively. The operator Wj appearing in the last line is obtained from the classical 
quantity by replacing the field functions with operators A afl (x), uo a {x), oJ a (x), ip aa (x) and ip a/3 (x) an d 
then taking the normal products with respect to |0). Each of the two proportion coefficients in (4.10) 
should make [0, 0, 0, 0, 0] equal to 1. It is evident that the time evolution of Hi is determined by i?o : 

Hi{t) = e iHot H/(0)e- iHot . (2.11) 

Therefore, the operator 

Texpj - i J d 4 xHi(x)} 

is the time evolution operator U(oo, — oo) in the interaction picture defined by Hq and the interaction 
/ d 3 xHi(0), namely 

Texpj - i J d i xH I (x)} = U(oo, -oo) , (2.12) 

and 



U(t, to) = e iff ot e -iH(t-to) e -iHoto ; (2.13) 

H = H a + J drxH^). (2.14) 
It follows that (2.10) can be written as 

4 R] [J,X,X,V,V] oc (0|C/( C o,0)(Tcxp{i|d 4 x[j^(x)A^(x)+x a (^a(^) 

+ uj h a (x) Xa (x) + Va(x)^(x) + ?!(x)%(x)] }) U(0, -oo)| ) , (2.15) 



The proportion coefficient should make [0, 0, 0, 0, 0] equal to 1. The time evolution of the operators 
A% , ■ ■ ■ are determined by H: 

wj(t) = e iSt uj a (0)e- iSt , 
uJ h a (t) = e iSt uj a (0)e- iSt , 
^(i) = e i5 Va(0)e- i5 *, 
^(i) = ^ a (0)e- & ■ 

Therefore H is the total Hamiltonian operator and A% (x), tjj(j:), ip%( x ) an d V'aO'O stand for the 

Heisenberg operators of the renormalized field functions. Thus (2.15) means that the renormalized Green 
functions can be expressed as the vacuum expectations of the time ordered products of the Heisenberg 
operators of the renormalized field functions. It should be emphasized that these field operators satisfy 
the usual equal time commutation or anticommutation rules. 



III. The scattering matrix 

One can choose to specify the initial condition at t = — oo or at t = +co for state vectors in the 
interaction picture and constructing the scattering theory. Therefore for an arbitrary state |*gb) in the 
QGB subspace there is a state \^' GB ) in this subspace so that 

C/(i,oo)|*' GS ) = C/(i,-oo)|* GB ). (3.1) 

Conversely, for an arbitrary state I^qb) £ QGB, there is a state \^gb) € QGB that satisfies this 
equation. Next, states U(t, ±oo)|*gb) are normalized just like \^gb) one thus has 

U(±00,t)U(t,±00)\* G B) = \*GB). (3.2) 

This and (3.1) also indicate that the operators U(±oo, Too) conserve invariant the QGB subspace. We 
therefore express the elements of the scattering matrix as 

1 £-o (0|f/(oo,-oo)|0) ' 1 ' 

where \i) and |/) are within the QGB subspace and are the eigen states of H . 
From (3.1) and (3.2), for an arbitrary state |*gb) G QGB one has 

U{0,-oo)\* GB ) = U(0,oo)\9' GB ) , 
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£/(oo,-oo)|* GB ) = |*' GB ) , 
C/(-oo,oo)|f GS ) = |* GB ). 

Thus, 

U(-oo,oo)t7(oo,-oo)|* GB ) = f/(-oo,oo)|tf GB ) = |* GB } . 
Similarly for arbitrary I^qq) E QGB, 

C/(oo,-oo)C/(-oo,oo)|* GS > = |*' GB ). 

These ensure the unitarity of the scattering matrix. 

With the help of the time evolution operator in the interaction picture we can also define the operators 
of so-called asymptotic fields which will be denoted by A aftt i n (x), ^4 a ^,out(a;), ipa, in (x) and ip a ,out{x), 
These are free fields, and when —t is very large each in-operator proportional to the Heisenberg operator 
of the renormalized field function. Actually each proportion coefficient should be equal to 1 in order not 
to destroy the usual equal time commutation or anticommutation rules. On the other hand when t is 
very large each out-operator is equal to the Heisenberg operator of the renormalized field function. Using 

A h a ^(x) = U(0,t)A^(x)U(t,0), 
rt(x) = U(0,t)<p a (x)U(t,0), 
^(x)=U(0,t)uj a (x)U(t,0), 

and noticing that when \t\ is large enough, U(0,t) and U(t,0) are independent of t, one gets 

A a ^ in {x) = U(0, -oo)A a ^x)U(-oc, 0) , (3.4) 
lk,in(a;) = U(0, -oo)ip a (x)U{-oo, 0) , (3.5) 

A-afx, out 

) = U(0,oo)A a „(x)U{oc,0), (3.6) 
V'a.outW - U(0, oo)^ a {x)U(oo, 0) . (3.7) 

Since the asymptotic fields are free ones, these are valid for any t. The formulae for the other fields are 
similar. 

If the in-states and out-states are defined by 

f/(0,-oo)|0) 



|0,in) 



|0,out) = 



v/(0|(7(oo,-co)|0) ' 
f/(0,oo)|0) 



v/(0|t/(oo, -oo)|0) 



^(0,-00)10 



\i, out } 



V<0|£7(oo,-oo)|0> 
U(0,oo)\i) 



^/(Ol^oo.-ooJlO) ' 
then (3.9) can be written as 



S/i = lim(/,out|i,in) = lim(/,in|f7(0,-oo)Z7(oo,0)|i,in) . (3.8) 



IV. An expression for the scattering matrix elements in terms 
of the truncated renormalized Green functions 

In the following we will derive a formula to express the elements of the scattering matrix in terms of 
the truncated renormalized Green functions. Define the following operator 

U{J,X,X,V,V) = T (jj(oo, — oo)cxp|i J d 4 x[j^(x)A a ^(x) + Xa( x ) LJ a(x) 

+uJ a (x) X a(x) +rJ a (x)Mx)+^a{x)Va(x)]}^) . (4.1) 

We thus have 

17(0,0,0,0,0) = [/(oo,-oo), 

and 

{0\U(J,x,X,%v)\0)/((0\U(oo,-oo)\0))=Z^ [J, x, X , %v}- (4.2) 
Analogous to (2.10), we also have 

u(J,x,x,v,v) 

= exp(-i f d 4 xH I ( 5 v . J*, w J x , ^ -J, v 7 .<f x , n ) 1 U W (J,% X,V,v) , (4-3) 
[ J \iSJ(x) iS x (x) (-i)Sx(x) i6r)(x) {-i)5r](x) ) J 

where 

U W (J,X,X,V,V) = Tcxp[iJ d i x[j^(x)A a ^x)+x a (x)u; a (x) 

+uj a {x) X a{x) +r] a {x)tp a {x) + tp a {x)r] a (x)] } . 

Define Hi(x) and C' GF {x) so that 

/off = - j d^Hjix) + j d 4 x{£ GF {x) + £ {x)} , 
C' GF = Cgf + \(d»A a ,) 2 = (d»A a ,) 2 . 
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We can write (4.3) as 
U(J,x,x,% v) 

= exp(-i f d A xUi( 5 r , 7 x , .c-/ T ' 7 \ ) 1 ( J ' X. X, ^ ??) . ( 4 - 4 ) 

t J \iSJ(x) iS x (x) (-i)5x(x) iSr](x) {-i)Sr](x)J J 

where 

P (0) WX,X,'),1)= exp|i| d 4 x£' GF (r^^}uW(J,x,x,V,v)- (4.5) 

Noticed that (0|£/(°)(J,x,X>'7> J ?)|0) ^ proportional to the generating functional for the Green func- 
tions defined by the effective Lagrangian £ + £gf and the latter can be explicitly expressed as 

zf\j,x,X,%V\ = ex p{-^/ J d 4 yJZ(x)Df v (x-y)JZ(y)} 

xexp|-i J d A x J d i yx a {x)C ab (x - y) X b(y)} 

where iD^,(x — y), iC a fc(x — y) and iS , a b(a; — y) are the propagators. Particularly, one has 

°?M = fc2 _M 2 + ie {^ + K - ^ + ' (4 ' 6) 

[f(fc)" 1 ]^ = { - (fc 2 - M 2 + ie)«T + (1 - |)fc"fc"}5 6 • (4-7) 
As for the operator U^°\j,x,x,r}y r ))^ one has 

u^(j,x, x ,v,v) 



-11 v v 



oc:exp{-y d A x J d A y[A a ^x)[D(x-y) iab SJ ^ y} 

6 S 



■+<%(*) - yrXa^T + [C(x y)-\^ a {y)-^—) } : Z^[J,X,X,V,v] ■ 



Substituting this in (4.4) we get 

C/(J,X,X,??^)/((0|C/(oo,-oo)|0)) 



= :exp{- J d A x j d A y(A a ^x)[D{x-y)- 1 \ 



ab 6JZ(y) 



s s 
+^(x)[s(x v )-\ , OQ ^y + [Si' yr% , a JMs^) 

+ zj b (x)[c(x yr 1 ]^^ + [Ci* y)-\Mv)^)) } : 4 R] ^x,x,%v] ■ (4.8) 



Denoting by iD^ u (x - y), \C ab {x - y) and iS ab {x - y) the full propagators denned by [J, \, X, V, v], 
and changing the sources in (4.8) into 

J?{x) = J d*x> J a ,(x') i [D(x'-x)-r a D , 

V'aa( x ) = - J d*x' T] a , a ,(x')i\S(x' -X)- 1 ]^^^ , 
V'ac( x ) = - J dVi[5(x-x') _1 ]aa, a 'a' Va'a'(x'), 
Xa( x ) = - J d'x'x'ai^ilCix'-x)- 1 }^, 

X' a {x) = - J d A x' i \C{x-x')-% a , Va ,(x') , 
we can find, in the limit as £ — > 0, that 

s,i = (/l : exp { /" 4x ( " ^ (8) ii) +f « (j) sra " ^ (a) (-i)«L(*) 

+IJa(x) iSx~(x) ~ (-i)Sxa(x) ) } 1 \ 1 ) X[R] ^X>X,V,V] \j =Y=x =rj^=o . ( 4 - 9 ) 

where [J, x, ^, 17] stands for the generating functional for truncated renormalized Green functions 
[10]. This is the formular we need. In the above manipulation, we have taken into account the fact that 
the momenta of the particles in the initial or final states are on the mass shell and for such momenta the 
renormalized full propagators become the free propagators. 



V. Concluding Remarks 

Starting from the path integral of the generating functional for the renormalized Green functions 
with the renormalized field functions as basical variables and using the unperturbed Hamiltonian operator 
that, under the Lorentz condition, describes the free particles of the initial and final states in scattering 
processes, we have formed a satisfying operator description and found the total Hamiltonian which 
determine the time evolution of the Heisenberg operators of the renormalized field functions. With the 
help of the time evolution operator U(oo, —00) in the interaction picture and of the Quasi-Gupta-Blculcr 
subspace formed by the physical initial-final states, we have been able to easily explain the definition and 
the unitarity of the scattering matrix and found a formular to express the matrix elements in terms of 
the truncated renormalized Green functions. 

We have also briefly described the asymptotic field viewpoint based on our operator description. It 
would be interesting to compare such a viewpoint with that of the traditional asymptotic theory [11 — 14]. 



10 



ACKNOWLEDGMENTS 



We are grateful to Professor Yang Li-ming for helpful discussions. This work was supported in part 
by National Natural Science Foundation of China (19875002) and supported in part by Doctoral Program 
Foundation of the Institution of Higher Education of China. 

Appendix The Gupta-Bleuler Subspace of Initial-Final States 

As was pointed out in section 2, the particles in the initial or final states are described by the following 
Lagrangian and Lorentz condition: 

A, = -\fM v FM»» + \M 2 A a ,A» + £*- (d^ a (x))d^ a (x) , (Al) 
d»A ati = , (A2) 

According to the operator description in the text of this paper we should keep the ghost term in equation 
(Al) and form the initial-final state subspace by following the Gupta-Bleuler quantization method. To 
this end we first use the Lorentz and modify (Al) into 

Co(x) = C GB (x) + ^M 2 A a) ,(x)AZ(x) + C+(x) - {d^ a (x))d»Lu a (x) , (A3) 

where the Gupta-Bleuler Lagrangian Cqb{x) can be written as 

C GB (x) = -^A av )(d^A:). 

Then we disregard the Lorentz condition and perform a canonical quantization with Co(x) and define the 
corresponding Hamiltonian operator Hq and field operators. Moreover we have to choose the so-called 
physical states and define a physical initial-final state subspace with the help of the Lorentz condition 
and the residual symmetry. We certainly expect the ghost particles to be excluded from appearing in the 
initial-final states. 

Assume that the operators of u) a (x) are anti-hcrmitian and therefore the operators of oJ a (x) is hermi- 
tian. Let w£\x), u>^\x) stand for uj a (x), iuj a (x) and Ili 1 ^, nl 2 \ n ajU stand for the canonical momenta 
conjugate to L0a\x), u^\x), A^ respectively. Thus the operators of these quantities in the Schrodingcr 
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picture are 



where 



n aM (x 



a£>(x 



ni 1} (x 
-i 2) (x 

ni 2) ( 



d 3 fc 



a aAt (k)e* x + at (k)e 



-ik-x 



VW V 2 

d 3 fc 



a aAl (k)e* k ' x - <(k)e- k - x 



V2|k|(27r) 3 
H)d 3 fc /|k| 



C a (k)e ik - x + ci(k)e- 4k - x 



vW V 2 

H)d 3 fc 



C (k)e lkx + Cl(k)e- lk - : 



d 3 fc 



C a (k)e 4kx -C](k)e 



(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 



H = -tr J d 3 fcft(|k|)<4(k)a ai ,(k) + J d 3 fc|k|{ci(k)a(k) + Cl(k)C a (k)}+i/°, (A10) 



n(|k|) = VM 2 + |k|2, 

a 0/1 (k), aL( k ')l = -guJab6 3 (k - k') , 



av(k),a bl ,(k') = <4 (k), o^(k') 



= 0, 



C„(k),C^(k /> 



C„(k),Cj(k / ) = <5 ab <5 3 (k-k') : 



C a (k),C b (k')J + = [C a (k),C 6 (k')J + = [C„(k),C b (k') 
C (k),C 6 (k')l =fc o (k),Cj(k / ) 



0. 



The restriction of the Lorentz condition on a physical state \^ p h) can be expressed as 

fe^a OM (k)|* pfc ) = (fc° = fi(|k|)). 
For a single particle state with the polarization vector (k) 

^(k,A)at M (k)|0> 

Eq. (16) gives 

V M (k,A) = 0. 



(All) 
(A12) 
(A13) 
(A14) 

(A15) 
(A16) 
(A17) 
(A18) 



According to this condition two transversal polarization states (A = 1, 2) and a longitudinal polarization 
state (A = 3) are allowed to be present for each k. 
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Since d fJ, d IJl ui a (x) = the theory we are treating is invariant under the infinitesimal transformation 

SA afl (x) = 5(dnuj a (x), Su a (x)=0, Scu a (x) = 0, 6ip{x) = 0, (A19) 

where 5( is an infinitesimal fermionic real constant. Similarly, since d^d^Uaix) = the theory is also 
invariant under the transformation 

6A a „(x) = 6(d^ a (x) , Su a (x) = , 5ZJ a (x) = , S^(x) = . (A20) 

This is the residual symmetry mentioned above. Under such a kind of transformations a physical state 
should become a equivalent state. Particularly e M (k, A)(5aj, M (k)|0) should be equivalent to zero, where 
<5aj lAt (k) is determined by 

5A ait (y) = 5C[d ll u a ] t=0 , (A21) 

or by 

M o/1 (x)=JC[^a] t=0 - (A22) 

From (A21), one has 

Sa^ (k) - i5(k» VWW\Cl (k) . ( A23) 

where the components of k^ are 

fc = |k| , h = ki . (A24) 

Similarly from (A22), one has 

<Jat M (k) = i5(k^^n/\k\cl(k) . (A25) 

It follows that C^(k)|0) and cj,(k)|0) should be equivalent to zero. 

In conclusion, the initial-final state subspace is formed by all the physical states. A physical state 
can only contain material particles, transversal polarization or longitudinal polarization gauge Bosons 
and not ghost particles. 
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